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We discuss models to calculate one- and two-neutron capture reac- 
tions on light nuclei. These are applied to calculate the reaction rates 
of 15 N(n, 7 ) 16 N, 16 N(n, 7 ) 17 N and 4 He(2n, 7 ) 6 He. The possible astro- 
physical importance is discussed. 



1 Introduction 

Neutron capture reactions on light nuclei play a role in various 
astrophysical scenarios. In the framework of Inhomogeneous Big 
Bang Models a high neutron flux can bridge the mass 5 and mass 8 
gaps. Subsequent neutron capture reactions may trigger a primor- 
dial r-process 0. Another site for neutron capture reactions is the 
high-entropy bubble formed during a type II supernova i'iS. Due 
to the photodisintegration at very high temperatures an a-rich en- 
vironment is created. When the temperature has dropped heavier 
elements can be built up mainly by a- and neutron-capture reac- 
tions. Again a critical question is how the mass 5 and mass 8 gaps 
can be bridged. 

In Section 2 we describe our model to calculate one-neutron 
capture reactions. This model is applied to the reactions 
15 N(n,7) 16 N and 16 N(n,7) 17 N. In Section 3 we will describe the 
theory of calculating a two-neutron capture reaction in a three- 
body model. We calculate the reaction rate of 4 He(2n,7) 6 He and 
compare the result with other works. In the final chapter we discuss 
and summarize our results. 
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2 One-Neutron Capture Reactions 

One-nucleon capture reactions on light nuclei are usually domi- 
nated by direct capture (DC) to bound states and resonant cap- 
ture to single, isolated resonances above the threshold. The level 
density above the threshold is normally low, i.e. between 0-10 lev- 
els per MeV. We calculate the cross section and reaction rate of 
the capture reaction with a hybrid model. The DC cross section 
is evaluated in a potential model. For the resonances we use the 
Breit-Wigner formula. We will discuss both contributions sepa- 
rately. 



2.1 Direct Capture 

The potential model is described by different authorsSi'B. We use 
real folding potentials as optical potentials!!!. 

The DC cross section of a transition to a bound state is deter- 
mined by the overlap of the scattering wave function, the bound- 
state wave function and the electromagnetic transition-operator. 
In most cases only El-transitions need to be taken into account. 

The total nonresonant cross section is determined by the 
direct capture transitions af c to all bound states multiplied with 
the single particle spectroscopic factors C 2 Si. 



^ = £(^K DC • (i) 

i 

The spectroscopic factors can be determined experimentally 
from other reactions, e.g., the spectroscopic factor of a (11,7)- 
reaction can be obtained from the (d,p)-reaction. Alternatively 
the spectroscopic factors can also be determined from shell-model 
calculations. 

^From the direct capture cross section we can obtain the reac- 
tion rate by integrating over a Maxwell-Boltzmann velocity distri- 
bution. We parametrize the non-resonant contribution to the rate 
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as 

N A < av > nr = A + BT 9 — CT° crn'mole'V 1 , (2) 

where Tg denotes the temperature in GK. The constant term repre- 
sents the s-wave capture contribution and the term proportional to 
T 9 the p-wave capture. All other contributions, i.e. higher partial 
waves and deviations from the conventional behaviours, are fitted 
in the third term. 



2.2 Single, Isolated Resonances 

The cross section of a single isolated resonance in neutron capture 
M processes s well described by the Breit-Wigner formula EM: 



where J and j t are the spins of the resonance level and the tar- 
get nucleus, respectively, E x is the resonance energy. The partial 
widths of the entrance and exit channels are T n and T 7 , respec- 
tively. The total width r tot is the sum over the partial widths 
of all channels. The neutron partial width T n can be expressed 
in terms of the single-particle spectroscopic factor C 2 S and the 
single-particle width r s p . of the resonance state EE 

r n = (C 2 S)T s . p . . (4) 

The single-particle width r s p . can be calculated from the scattering 
phase shifts of a scattering potential with the potential depth being 
determined by matching the resonance energy. 

The gamma partial widths T 7 are calculated from the electro- 
magnetic reduced transition probabilities B(Jj — > J/;L) which carry 
the nuclear structure information of the resonance states and the 
final bound states 0. The reduced transition rates are computed 
within the framework of the shell-model. 
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The resonant reaction rate for an isolated narrow resonance can 
be expressed a s00 



N A <av> T = 1.54 x 10 5 /i" 3/2 T 9 ~ 3/2 

x ^(c^exp (-11.605^/Tg) cnftnole^s -1 , (5) 

i 

where the resonance strength wy is defined as 

2J + 1 r n r 7 
W7= 2(^Ti)T— (6) 

and has to be inserted in eV in Eq. |5|. The resonance strength 
can be determined experimentally or derived from the calculated 
partial widths. 

2.3 One-Neutron Capture on N-Isotopes 

We will now consider one-neutron capture reactions on neutron- 
rich N-isotopes. We start with the reaction 15 N(n,7) 16 N which is 
known experimentally at stellar energies^. Therefore this reaction 
can be used as benchmark test for the validity of our model. 

The considered transitions for the direct capture of this reac- 
tion are listed in Table |l[ The four lowest states of 16 N can be 
described by a coupling of the l/2~ ground state of 15 N with a 
neutron from the ld.5/2- (resulting in the 2~ ground state and a 
3~ excited state) and the 2si/2-subshell (resulting in the low-lying 
excited states 0~ and 1~). In a shell-model description these states 
are good one-particle states. Therefore the spectroscopic factors 
should be close to unity. In fact the shell-model calculations re- 
sult in spectroscopic factors of about 0.9 for all statesB. However, 
the spectroscopic factors from a (d,p)-experiment are considerably 
lower by an average factor of around 1.8 0. The reason for this 
discrepancy is not known. In our DC calculations we use the ex- 
perimental spectroscopic factors. 
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Table 1: Considered transitions for the direct capture reactions on N-isotopes. 
Transitions with very small contributions are not included in the table. The 
Q-values are in MeV. 



reaction 


Q-value 




E x (MeV) 


transition 


C Z S 


ib N(n, 7 ) ib N 


2.491 




0.000 


P- 




0.550 






o- 


0.120 


P- 


- > 2s 1/2 


0.460 






3~ 


0.298 


P- 


+ld 5 /2 


0.540 






1" 


0.397 


P" 


^2s 1/2 


0.520 


16 N(n, 7 ) 17 N 


5.883 


1/2- 


0.000 


P- 


^ld 5 / 2 


0.589 






5/2~ 


1.907 


P- 


*ld 5 / 2 


0.207 






7/2~ 


3.129 


P- 


*ld 5 /2 


1.457 






5/2- 


4.415 


P" 


- > 2s 1/2 


0.921 



All transitions listed in Table [T] result from an incoming p-wave. 
The s-wave contribution which dominates at very low energies can 
be obtained directly by extrapolating the thermal absorption cross 
section with the l/u-law. 

The resonance parameters are listed in Table |2|. The total width 
of the 862 keV resonance of T = 15keV@ corresponds to the neu- 
tron width of the state. The 7 width was estimated with the rec- 
ommended upper limit as Y 1 = 4.2 eV@. ^From the shell model 
calculation we obtain a width T 7 = 0.455 eV which reduces the res- 
onance strength by about one order of magnitude. However, this 
resonance only has a small influence on the reaction rate at very 
high temperatures. Therefore this change of the resonance strength 
barely changes the reaction rate. 

The parameters of the direct capture reaction rate (see Eq. 0) 
are listed in Table |^. The resulting reaction rate agrees very well 
with the rate of Meissner et The rate is clearly dominated by 
the direct capture contribution. Both rates show good agreement 
with the experimental data. For a discussion of the experimental 
data we refer the reader to the paper of MeissnerB. 

The bound levels of 17 N are known from experiment. We used 
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Table 2: Adopted values for the resonance parameters for capture reactions on 
nitrogen isotopes. 



reaction 


E x 


.r 


E 


r n 


r 7 








[MeV] 




[MeV] 


[eV] 


[eV] 




[eV] 


iiJ N(n, 7 ) ib N 


3.360 


i+ 


0.862 


15 000 


0.455 




0.341 


16 N(n, 7 ) 17 N 


5.904 


7/2- 


0.021 


0.032 


4.80-10" 


-2 


0.015 




6.121 


5/2+ 


0.238 


1.2 


4.80-10" 


-2 


0.027 




6.325 


3/2+ 


0.442 


20 


5.46-10" 


-2 


0.022 




6.372 


7/2+ 


0.489 


20 


1.52-10" 


-2 


0.012 




6.373 


5/2+ 


0.490 


600 


0.110 




0.066 




6.470 


1/2+ 


0.587 


1750 


2.510 




0.501 




6.685 


3/2~ 


0.802 


12 500 


5.660 




2.263 




6.737 


7/2+ 


0.854 


70 


4.17-10" 


-2 


0.033 




6.835 


3/2+ 


0.952 


360 


0.478 




0.191 



Table 3: Parameters for the direct-capture contribution to the reaction rate. 





A 


B 


C 


D 


15 N(n, 7 ) 16 N 


3.18 


3783.4 


335.2 


1.716 


16 N(n, 7 ) 17 N 




3649.9 


437.5 


1.633 



spectroscopic factors calculated in the shell-model for the DC cal- 
culation. The important DC transitions are listed in Table [I|. Sev- 
eral levels are known above the threshold but without spin/parity 
assignments. It is impossible to assign the levels with the help of 
our shell model calculation. Therefore we use the shell-model en- 
ergies for the resonances. The adopted resonance parameters are 
shown in Table 0. The two resonances at 21keV and 238 keV are 
the main contributions to the reaction rate. 

Since there is no s-wave transition the parameter A of the re- 
action rate vanishes (see Table [3]) . We find a considerable enhance- 
ment of our rate compared to the calculation of Rauscher et al. Bin 
that work the resonance at 197 keV was given a hypothetical 5/2 + 
assignment. Since no resonance at lower energies was included, 
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our 21 keV resonance causes a strong increase of the rate at low 
temperatures. 

3 Two— Neutron Capture Reactions 

In this section we will discuss the possibilities to calculate a three- 
body reaction rate. We will apply the theory to the reaction 
4 He(2n, 7 ) 6 He. 

The reaction rate of 4 He(2n, 7) 6 He can be calculated as a se- 
quential two-step process. We call this method the 6 He = a+n+n 
approach. Here it is assumed that in a first step the unstable nu- 
cleus 5 He is formed via the reaction 4 He + n — > 5 He (negative 
Q-value Qi). The second step 5 He + n — > 6 He + 7 (positive Q- 
value Q2) is treated as a two-body problem in the initial and as a 
three-body problem in the final state. The values for spin-parities, 
Q-values and widths used in our calculations are given in Table |j. 
The total Q-value Q12 is the sum of the Q-values Qi and Q 2 - 



Table 4: Spin/parities, Q-values and widths (in MeV) used in our a+n+n 
calculations of the reaction rate for the reaction 4 Hc(2n, 7) 6 He. 



J 7r ( 4 He) 


J^He) 


J^He) 


Qi Qi Q12 Tc 


0+ 


3/2- 


0+ 


-0.89 1.87 0.98 0.76 



In this two-step model the reaction rate of 4 He(2n,7) 6 He is 
given by 



Nl c (2n 4 He) = 2 / / dE 1 dE 2 



h d < n 4 He > d < n 5 He > 



r(^i 



dE 1 



dE, 



(7) 

The energies E\ and E 2 denote the center-of-mass collision energies 
of the first (second) neutron with the 4 He ( 5 He) nucleus. The 
width of the intermediate state, T(Ei), is energy-dependent. The 
differential rates of the two steps represent the product of the cross 
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section with the Maxwell-Boltzmann distribution: 

d<av> ( 8 \ f 1 \ 3/2 /rnXrn / E 

-) cr(E)Eexp(- 



dE V^/V \kT J \ kT 

The whole procedure of calculating the cross section is described 
in detail by Efros et a/0 and by Balog hi. 

The numerical results for the reaction rate (2n 4 He) calcu- 
lated in this model are fitted to the expression 

N% (2n 4 He) = £ cm 6 s" 1 mole -2 (9) 



in the temperature regions 0.1 < Tg < 1 and 1 < Tg < 3. ^From 
a comparison with the situation of small resonances it should be 
possible to perform a fit to a function ~ aiTg 2 exp(— a 3 /T 9 ), but it 
turned out that such a fit leads only for Tg > 1 to useful results. 
Therefore, we decided to use the parametrization of Eq. ([]), which 
reproduces the calculated data with an error less than 0.25% in the 
whole temperature range. The parameters are given in Table [5[ 

In the work of Danilin et al El the wave functions for ground 
and scattering states of the halo nucleus 6 He are calculated in an 
a+n+n three-body model and are used to predict the strengths 
of nuclear and electric dipole excitations. From the astrophysical 
point of view we are especially interested in the electric dipole 
strength function dB E1 /dE 7 , which is directly connected to the 
cross section for photodisintegration a^ 1 of the nucleus 6 He into 
4 He and two neutrons by0 

°i ~ —hc-^ (10) 

where E 1 = E + S^n is the photon energy, which can be written 
as the sum of the two-neutron energy E = E 1 + E 2 and the two- 
neutron separation energy, S2 n = Qn- We refer the reader to the 
work of Danilin et a/@ about details of the calculation. 
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Table 5: Fit-parameters a, for the reaction rate N% (2n 4 He) calculated in our 
6 He = a+n+n approach and from the El-strength obtained in the work of 
Danilin. 







a + n 


+ n 






Danilin et al 






0.1 < T 9 < 1 


1 < T g < 3 




0.1 < T 9 < 1 


1 < Tq < 3 




a 


2.75605E 


- 16 


-4.31090E - 


11 


-4.94120E 


- 11 


2.27413E - 


10 


Ol 


2.88001E 


- 12 


1.95298E - 


10 


2.86628E 


- 09 


-2.40420E - 


10 


«2 


4.41655E 


- 12 


-3.29920E - 


10 


-4.90850E 


-08 


5.64850E - 


09 


«3 


-1.35270E 


- 12 


2.77932E - 


10 


3.46189E 


- 07 


-2.23220E - 


09 


(I4 


1.68846E 


- 11 


-1.02300E - 


10 


-1.08100E 


-06 


3.69620E - 


10 


a 5 


-2.45810E 


- 11 


1.80468E - 


11 


1.87337E 


- 06 


-2.62480E - 


11 


a 6 


1.41336E 


- 11 


-1.26680E - 


12 


-1.85880E 


-06 


4.500240E - 


13 


a 7 


2.08214E 


- 11 






9.89903E 


-07 






as 


-2.67550E 


- 11 






-2.19670E 


-07 






ag 


8.22640E 


- 12 















In Fig. H the reaction rates N\ (2n 4 He) for the formation of 
6 He calculated in the ct+n+n approach and deduced from the El- 
strength function of Danilin et a/HBEJ are depicted. Furthermore, 
the results obtained in previous studies of Fowler et afl and Gorres 
et a/il are shown (we have used the constructive rate of that pub- 
lication). Our results obtained in the a+n+n approach show a 
relatively good agreement with those obtained by Gorres il, who 
also used a two-step model with a simple direct capture calculation 
for the second step. Below Tg = 0.8 no data for the rate N\ (2n 4 He) 
are available in the work of Gorres et a/H. Our calculations are 
larger by more than three orders of magnitude compared to the cal- 
culations of Fowler et a/il. This is due to the fact that in the work 
of Fowler it was assumed that the second step ( 5 He + ra — > 6 He + 7) 
proceeds via a resonant state in 6 He near the threshold followed by 
the emission of E2-radiation, which is much more less likely to take 
place than the nonresonant El direct capture of a neutron into the 
ground state of 6 He. 
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lO" 13 

10 -14 _J 

lO" 15 

10 -16 



a + n + n-approach 



Danilin et al. 

Fowler et al. 

• gorres et al. 



T 1 1 1 1 — I — I — I I I I I I I I I I I I I I I I I I 



0.1 



1 



n 



Figure 1: Comparision of different rates (2n 4 He) for the formation of ' 



The most remarkable fact is that the rate deduced from the 
El-strength function of Danilin et a/Sill is for T 9 > 0.2 approxi- 
mately two orders of magnitude larger than our results obtained in 
the a+n+n approach. A possible explanation of this enhancement 
is that in our a+n+n calculations only the sequential process is 
taken into account whereas in the calculations of Danilin et al III 
the simultaneous decay ( 6 He + 7 — > 4 +n+n) is also included. Fur- 
thermore, it is not possible to describe the threshold behaviour of 
the cross section a^ 1 (with the help of the reciprocity theorem) 
very well for broad resonances in a two-step model. 

4 Summary and Discussion 

We calculated reaction rates for one-neutron capture on 15 N and 
16 N and two-neutron capture on 4 He. The rate for 15 N(n,7) 16 N is 
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in good agreement with both previous calculations and experimen- 
tal data. In general the neutron capture rates on stable targets 
are quite well known. The reaction 16 N(n,7) 17 N is an example of 
a capture reaction on an unstable target. We find a considerable 
enhancement to a previous calculation. This enhancement can be 
also be observed at various other capture reactions on unstable 
targets in this mass range 0. This fact could influence the reac- 
tion path both in the nucleosynthesis of Inhomogeneous Big Bang 
Models and in the alpha-rich freeze-out of type II supernovae. 

For the reaction 4 He(2n,7) 6 He we also find a strong enhance- 
ment of the rate compared to the calculation of FowlerE3. But from 
the calculation of the inverse photodisintegration rate by Danilin 
et alE we deduce an even much higher rate. This is probably due 
to the fact that this rate also includes the simultaneous capture 
of two neutrons. But even with this enhanced rate the path via 
4 He(2n,7) 6 He is dominated by the reaction 4 He(an,7) 9 Be at con- 
ditions typical for the alpha-rich freeze-out. This is due to the 
effective destruction of 6 He through photodisintegration. 
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